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The many electron correlated scattering (MECS) approach to quantum electronic transport was 
investigated in the linear response regime [I. Baldea and H. Koppel, Phys. Rev. B. 78, 115315 
(2008)]. The authors suggest, based on numerical calculations, that the manner in which the method 
imposes boundary conditions is unable to reproduce the well-known phenomena of conductance 
quantization. We introduce an analytical model and demonstrate that conductance quantization is 
correctly obtained using open system boundary conditions within the MECS approach. 
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I. INTRODUCTION 

Electron transport calculations in the linear response regime for a method referred to as the stationary Wigner 
function (SWF) method have been reported [I. Baldea and H. Koppel, Phys. Rev. B. 78, 115315 (2008)], as a 
shorthand we will refer to their presentation as BK. In fact though, their paper is a comment on the validity of 
a method introduced by Delaney and Greer [T] for treating open system boundary conditions for correlated many- 
electron problems, or many-electron correlated scattering (MECS) [HO [3]. For reasons that will hopefully become 
clear, we refer to the method under discussion as MECS [J. BK express several criticisms of the MECS method, but 
their primary objection and conclusion is 

"... the manner in which it" (open system boundary conditions through the Wigner function) "was imposed in" 
[P. Delaney and J.C. Greer, Phys. Rev. Lett. 93, 036805 (2004)] "turns out to be inappropriate. It misses the 



'Electronic address: jim.greer@tyndall.ic 



^Electronic address: p. dclaneyi9qub.ac.uk 
•^Electronic address: jgeorgios.t'agas@tynd all.ie| 



2 

fact that, in accord with our physical understanding, the current flow is due to an asymmetric injection of electrons 
from reservoirs into the device, and that injected electrons are very well described by Fermi distributions with different 
chemical potentials. Moreover, as results from the analysis at the end of Sec. VI" (of the BK paper), "unfortunately 
there is no simple remedy of the SWF method; the modification of the boundary conditions in the spirit of" Delaney 
and Greer "such as to account for a nonvanishing chemical potential shift does not yield the desired improvement. " 

The conclusion in BK is reached after numerical calculations on a test system designed to model independent 
electron and correlated electron transport across a quantum dot. As we will demonstate, their conclusion incorrectly 
presumes an asymmetric injection of momentum is necessary to describe reservoirs or electrodes in quasi-equilibrium, 
and likewise incorrectly assumes there is no chemical potential difference when applying the open system boundary 
conditions with the MECS procedure. We will show that the Wigner boundary conditions as expressed for MECS 
calculations, or equally in previous works [5], is consistent with the application of a 11 non-vanishing chemical potential 
shift" . The authors in their paper appear to be confusing the application of open boundary conditions when using 
either energy or momentum distributions, and we will clarify that a relative shift in energy to the reservoir energy 
distributions does not imply a shift for the corresponding momentum distributions emitted into a scattering region 
for simple models of electrode behaviour. 

In sect. [TTJ we give a brief overview of the MECS method. This is followed in sect. |III| by introduction of a resolution 
to the apparent conundrum expressed by BK: we analyse momentum distributions for electron reservoirs or electrodes 
represented by Fermi-Dirac distributions and observe that the momenta distributions do not change with applied 
voltage bias, or equivalcntly with a chemical potential imbalance applied between electrodes. Asymmetric injection 
implied by BK's conclusion in relation to the momentum distributions from the electrodes is not consistent with 
a MECS or a Landauer description of electron transport. In sect. |IV| a calculation of conductance quantization 
for a system with the electrodes represented by parabolic energy bands in one spatial dimension is given using 
the MECS construction. To clearly identify how the boundary conditions can be treated in this case, electron 
interactions are neglected and the method reduces to a many-electron scattering problem for non-interacting electrons. 
This approximation has the advantage of highlighting the simplicity of the MECS formulation as well as clearly 
demonstrating its consistency with standard formulations of electron scattering, while avoiding issues related to 
specific numerical approximations or questions related to specific electronic structure implementations. The boundary 
conditions as formulated in MECS applied to the model reproduces the well-known result of conductance quantization. 
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II. A BRIEF INTRODUCTION TO THE MECS APPROACH 

The proposal behind MECS as introduced in ref. [I] is to constrain a many-electron wavefunction on a scattering 
region, or specifically, the many-electron density matrix (DM) 

PN = |*JV >< *iv|, (1) 

in a manner satisfying open system boundary equations [5]. As open system boundary conditions are commonly 
expressed in the language of single particle theories, it is useful to consider the Wigner transform of the one-electron 
reduced density matrix (RDM) associated with the many-electron DM 

Pn — * pi, with Trpi = N 

r+oo 

fw{q,P) = / dr exp(-ipr/K) Pl (q-r/2;q + r/2), (2) 



where here the transform is written for one spatial dimension r and {p, q} are Wigner phase space variables. The MECS 
approach is the recognition that open system boundary conditions can be applied to correlated systems through the 
one-body RDM through use of the Wigner transform, allowing appropriate conditions to be enforced at the boundaries 
of a scattering region. In practice, the Wigner distribution function is used to constrain the momenta flow out of the 
electron reservoirs and into the scattering region. The momentum expection value can be written as 
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<p> ^2irhj d P d i pfw(q,p)- (3) 

Eq. [3] highlights the analogy of the Wigner quantum phase space distribution to a classical probability distribution 
function. However, unlike a classical probability distribution and as is well known [51 El El [5], the Wigner function 
is not everywhere positive as a consequence of the Heisenberg momentum-position uncertainty principle. However, 
in regions where fw behaves approximately classically, the Wigner phase representation allows us to assign meaning 
to phrases such as "electrons in the left or right reservoir" , "momentum of an electron emitted from a reservoir" , 
or "a reservoir is locally in equilibrium". Within this context, the net momentum flow out of the left electrode is 
approximated as 
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+00 



and similarly for the right electrode 



Pi = ^ J dpp fw(quP), (4) 



1 f° 

Pr= 2irhJ d PPfw{lr,P) (5) 
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where qi and q r are appropriately chosen. Clearly this approximation is dependent upon how well fw describes a 
classical probability distribution function. For metal electrodes, where electrons can be well approximated by the 
free electron model, this assumption is usually justified. It is also noteworthy that the Wigner reduced one-particle 
function, as a function of energy defined in the Wigner phase space, tends rapidly toward the Fermi-Dirac distribution 
with increasing number of particles in a confining potential [S] . 

In our calculations to date, three-dimensional electrodes are considered [TJ [3J []Jjl HI] . To simplify the analysis, fw is 
integrated over the in-plane coordinates within a cross section of the electrodes. The net momentum flow out of both 
electrodes is constrained to this equilibrium (V = 0) value [HH^]- To consider this procedure further, we examine 
models of electrode behaviour in quantum transport theories. 

III. FERMI-DIRAC RESERVOIR ENERGY AND MOMENTUM DISTRIBUTIONS UNDER APPLIED 

VOLTAGE AT ZERO TEMPERATURE 

In fig.JIJi), a pictorial representation of open system boundary conditions is shown for electrodes described by two 
parabolic energy bands (free-electrons with effective mass m* in one dimension and for a Fermi-Dirac distribution at 
temperature T = 0). The energy levels eu = in the left and right electrodes are filled to the Fermi energy ep ; 

similarly momentum states are filled to the Fermi momentum /cf. In fig.]!]}), a potential energy difference is introduced 
between the left and right electrodes shifting the bottom of the energy bands with respect to each other by an amount 

denoted eV. This results in a shift to the energies in the right electrode by > ^2m* e ^ ( a symmetric split in 

the voltage between electrodes does not alter our discussion). The shift of the right electrode energy states describes 
the chemical potential imbalance between the reservoirs and will necessarily be accompanied by a voltage drop, or 
equivalently an electric field across the scattering region. In fig. [TJ:), the momentum distributions corresponding to 
Fermi-Dirac energy distributions with and without applied voltage in fig. [TJi) and fig.jljD), respectively, are shown. The 
same momentum distributions are obtained with or without application of a voltage difference between the left and 
right electrodes. It also follows that similar considerations hold for the case of non-zero temperature. It is important to 
highlight at this point that the MECS proposal [TJ for treating quantum electronic transport is completely compatible 
with this picture of boundary conditions on the electrode regions. 

When working with the one-electron RDM obtained from a correlated iV-electron density matrix, it is not possible to 
unambiguously define single electron energies. It is therefore of advantage to constrain the total incoming momentum 
to model the action of the electrodes while determining the correlated many-electron density matrix on the scattering 
region. As the incoming momentum distributions are the same for electrodes in equilibrium or in local equilibrium, 
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constraining the net momentum flow from the electrodes is not sufficient to drive the electrodes away from equilibrium 
with respect to each other. It is standard practice in many electron transport methods, including non-equilibrium 
Green's function (NEGF) techniques, to introduce an external electric field to model the action of the electrodes on 
the scattering region. To understand the role of the external field in transport calculations, it is useful to consider the 
sketch of two metallic electrodes as presented in fig. [2] A simple band or independent particle model is a remarkably 
good approximation to the electronic structure of metals allowing our discussion of the parabolic bands to be extended 
directly to consideration of realistic models of metal electrodes. Within fig. [2] the fact that the two electrodes are 
not in equilibrium with respect to one another is denoted by the different left fi^ and right chemical potentials. 
The chemical potential imbalance introduces a difference in the charge density between the left and right electrodes. 
However, electrostatic screening is efficient in metals and for the quasi-equilibrium regions of the electrodes, the 
electric field is zero within the electrodes or equivalcntly, the voltage is constant within a short distance into the metal 
electrodes. Thus all of the voltage drop is across the scattering region plus the screening length into the electrodes. 
Typical screening lengths in metals are of the order of 0.1 nm. Hence any charge imbalance in the electrodes resides 
at the surface of the metal and the opposite polarity of the surface induced charges between the electrodes gives rise 
to an electric field across the region situated between the electrodes; a situation depicted in fig. [2] as field lines between 
the electrode surfaces. In most transport calculations, charges in the electrode and scattering regions are solved for 
self-consistently allowing a molecular tunnel junction to polarize in response to this external electric field. In MECS, 
charges re-arrange due to minimization of the energy with respect to the many-electron wavefunction subject to the 
open system boundary conditions and the external electric field. The voltage can then be extracted as the combined 
field arising from the applied field and polarized charge distribution in the scattering region; see for example |13) . 

The model of electrode behaviour we are describing is consistent with other formulations of quantum electronic 
transport. In this regard, it is worthwhile to mention the work of McLennan, Lee and Datta |14j and in particular 
their fig. 3. As pointed out by the authors, a change in the chemical potential in the electrodes is accompanied by a 
shift in the electrostatic voltage. In a metal, the case considered for MECS calculations to date, the chemical potential 
and electrostatic voltage changes are nearly identical and the effect of the applied voltage can be accounted for as a 
shift in the electrode bands with respect to each other. In the linear response regime, the shift due to the electrostatic 
voltage is sometimes neglected |14j but formally should be included as indicated schematically in our fig. [I] and in 
McLennan et al's fig. 3. 

The point which we would like to emphasize is that the presence of voltage drop on the scattering regions ensures 
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that the single electron energies in the electrodes are shifted by the applied voltage, but does not alter the momentum 
distributions emitted from the electrodes. The conclusion in BK that an asymmetric injection of electrons is needed 
to obtain a current is incorrect, if injection refers to incoming electron momentum distributions, i.e. 'current injected' 
from the electrodes. What is required is asymmetric scattering for injected electrons, and this is provided for by the 
asymmetric electric field profile or equivalently, the chemical potential imbalance generated across a molecular tunnel 
junction. 

IV. APPLICATION OF MECS TO A SINGLE PARTICLE MODEL 
A. Boundary conditions 

MECS was originally formulated for interacting many-electron systems with a Hamiltonian operator given as 

H N = T +V1 + V2 (6) 

where T is the sum of N one-electron kinetic energy operators, V\ is the sum of the external one-electron potentials 
and V2 is the sum of electron-electron interactions on the scattering region [T]. As mentioned, with the MECS 
approach an external electric field can be applied to model the chemical potential imbalance between the reservoirs 
and this term may be included into V\. To arrive at a single particle model to be used in our analysis, the electron- 
electron interactions are 'switched' off and, for simplicity, all other external potentials other than the chemical potential 
imbalance between the electrodes are also 'switched' off. As electron-electron interactions are not treated, the resulting 
many-electron system consists of N non-interacting electrons each described by a single electron Hamiltonian operator. 
Of course, there is no advantage to apply the MECS method to a non-interacting electron model, but we demonstrate 
in this case that the MECS method is consistent with the usual formulation of quantum mechanical scattering. 

We consider a scattering region of length [—L/2, L/2] in the absence of a potential and write plane wave eigenfunc- 
tions on the scattering region 

ip n (x) = -^=exp(ik n x) 

k n = 2nn/L ; n = ±1, ±2, ±3, . . . (7) 

Left going and right going states are filled to a Fermi momentum kp — npAk. In the absence of the application of a 
voltage, the current is given simply as 

I[V = 0} = -£^ E k n \=0, (8) 

\n(=l n r = l J 



where rt/ and n r denote states associated with the left and right electrodes, respectively. In one dimension, the electron 
current and current density are equivalent, and the factor A indicates our choice of normalization. The model consists 
of left and right propagating plane wave states with the net current summing to zero. In this representation, the 
density matrix is diagonal with n„ „< = 5 nn > for both n,nf < rip and n n n i = otherwise. This allows the density 
matrix to be constructed from the first Hf states incoming from the left and right to be written as 



p (x, x') = — exp [ik„(x - x')] + — ^ CX P MM^ - x')] . (9) 

n— 1 n— 1 

Introducing the Wigner transformation term by term, the resulting Wigner distribution function is readily found to 
be 

n F 2 ^ n F 2 ^ np 

n— 1 n— 1 n=l 

with S(p) the Dirac delta function. Strictly speaking, as we consider wavefunctions that are only non-zero on the 
scattering region, the delta functions should be replaced by sine functions of the form: ^j^^jj^ sin [(fc ± p/H)L/2] 
that approach delta functions for large L. In the case of large L, the model as described corresponds to fig. |TJt) . 

Next, a potential step is introduced at x = to drive the system out of equilibrium and allow for a net current 
flow. The exact form of the scattering potential is not critical to the following argument, but for ease of presentation 
we assume that the potential is varied over a small region I « L allowing us to approximate the difference between 
the left and right electrodes as a potential energy step. For this case, the solutions to the one-electron Schrbdinger 
equation may be written in scattering form 

ip n (x) — exp(ik n x) + rexp(—ik n x) x < 

4>n(x) = texp(ik' n x) x > 0. (11) 

For spatially varying potentials centered at x = satisfying I << L, the asymptotic wavefunctions will likewise satisfy 
the scattering form and the following development remains valid with minor modification. By implication, the energies 
for electrons entering the scattering region from the right are shifted by an amount given by the scattering potential 
height ^2m» ~^ ^2m* e ^ ' ■ However the incoming momenta, as previously discussed, are unchanged. We apply a 
voltage greater than the level spacings at the bottom of the reservoir conduction band where the energy density of 
states is greatest. With the introduction of the potential, the model for electrons entering from the left is the one- 
dimensional quantum mechanical scattering problem with a step-up potential. For electrons entering the scattering 
region from the right, the problem is for an electron incident on a step-down potential. This is shown schematically in 



fig. [TJd) and the solution is well known for both cases and may be expressed in terms of the transmission coefficients 
T. In this case, we can write the current for the system as 



I[V^Q] = —^—y d [fi(k n )Ti{k n ; V)k n - f r (k n )T r (k n ; V)k n ] , (12) 
m*L 

n=l 

with Ti,T r the transmission coefficients for electrons incoming from the left and right, respectively. Note that fi, r (k n < 
kp) = 1 and /; r {k n > kp) = for our example. Time reversal symmetry requires that the left and right transmission 
functions for a given single particle energy e are equal T;(e) = T r (e), but the energies for the left and right states of 
equal momentum are not equal in the presence of a voltage. This is seen by re-writing the transmission as functions 
of energy 



K, (13) 



resulting in a net current flow with application of voltage. 

It can be shown that the Wigner transform of the reduced density matrix constructed from the scattering wave 

functions satisfies the same open system boundary conditions as the zero voltage (plane-wave states) solution in the 

large L limit and as we will demonstrate, approximately for finite values of L typically used in numerical studies. The 

electron reservoirs in this case are the regions outside of the central scattering site. The density matrix with V =/= 

remains diagonal and we can again consider the Wigner transform term by term 

f + °° r r 

fn,nW>P) = / dr exp(-ipr/h)tp n (q+ -)%l)* n {q - -). (14) 



where now ip n are wavefunctions of the scattering form eq. 11 on [-L/2, L/2). We consider a point to the left of the 
scattering region qi = —L/A and find for the Wigner distribution function 

2lT 2lT o 47T 

fn,n{<luP) = -^ 5 (P~ hk n) + -j-r S(p + hk n ) + —r cos(2k n q)S(p) (15) 

where in eq. [15] r is the amplitude of the reflected component of the scattering wavefunction and the normalization is 
fixed to that of an incoming plane wave. The Wigner phase space density at qi consists of the incoming momentum 
term at p — +hk n , the reflected momentum term at p — —hk n , and a zero mode term p = as discussed in ref. [6]. 
The Wigner transform is calculated at ±L/4, in the middle of the electrodes, to avoid coupling to the regions outside 
of [—L/2, L/2] and to avoid interaction between the electrodes as voltage is applied. If the point where the Wigner 
function is to be constrained is too close to the boundaries ±L/2, the density matrix, as can be seen from the argument 



of the Wigner transform eq. 14 ip n (±L/2 + |)-0*(±L/2 — |), will be zero. The vanishing of the density matrix in 
this case is an artifact arising from truncating the wavefunctions outside of the scattering region. If the point where 
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the Wigner function is to be constrained is calculated too close to the scattering region, as voltage is applied the 
incident and reflected components of the scattering wavefunction mix with the transmitted component, or in other 
words the two electrodes couple. To avoid coupling the electrodes, the Wigner function should be calculated within 
each electrode, but in a region avoiding interaction between the electrodes as a voltage is applied. 



Again, as for eq. |10| for finite L the delta functions in eq. 15 should be replaced by sine functions of the form 
Tr(k±p/h) sm [(^ ^ For large L when the sine functions well approximate delta functions, eqn 
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and 



15 



satisfy the same condition for incoming momenta states. For finite L, the sine functions centered at ±fc, can overlap 
and the incoming momentum states can differ between the V — (plane-wave) and V ^ (scattering) states. In 
fig. |3^l) , we have calculated the Wigner function for a set of plane wave states incoming from the left and the right 
with a Fermi momentum kp chosen to correspond to that of gold electrodes and a scattering region of L = 2 nm. A 
step potential of V = 0, 1, 2 Volts is applied in the center of the scattering region and the resulting Wigner functions 
are displayed at qi = —L/A for each value of the applied voltage. These parameters have been chosen to compare our 
analysis to typical calculations for molecular electronics, and in particular the calculation presented in ref. pQ. As 
voltage is applied, the incoming electron distributions are not exactly equal to the V = distribution, but agreement 
is very close particularly for states near kp which provide the largest contributions to the current. The outflow of 
momentum is changed as a result of the scattering off the potential barrier and this difference between the equilibrium 
(V = 0) and non-equilibrium distributions (V =/= 0) reflects the difference in chemical potential between the electrodes. 
In fig. [3j>) , the length of the region is taken to be L — 20 nm and as seen the incoming momentum distributions, 
even with wavefunctions only defined on the region [—L/2,+L/2], remain essentially the same for V — 0,1,2 volts. 
The momentum distribution drops sharply at kp well approximating the distribution shown in fig. [lb). In fig. [3b ) 
and d), the Wigner distribution at q r = +L/4 for L = 2 and 20 nm, respectively, is given. In subsequent discussion, 
we assume a large value for L, but as fig. [3] indicates, the considerations apply well to electrode lengths as small 
as 1 nm. Also, the lowest rip states have been occupied corresponding to a T — distribution. However, there is 
nothing in the analysis that precludes the V = solution to be set to a thermal distribution and to constrain the 
solutions to the thermally occupied incoming states as a voltage is applied. If the scattering states with V ^ are 
constrained to satisfy the Wigner function determined from the V = wavefunctions, the correct solution to the one 
electron Hamiltonian in the presence of a potential on the scattering region will be obtained in the large L limit and 
approximately for smaller values of L. Indeed it is observed that in the single particle case, constraining the incoming 
momentum inflow via the Wigner function implies solving the one-electron Schrodinger equation for a specific value 
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of incoming momentum p — hk n , and otherwise results in the standard textbook presentation of one-dimensional 
quantum mechanical scattering. 



B. A transport calculation with the single particle model 

We again consider introduction of a small potential step to shift the energies of the states incoming from the right, 
as depicted in fig.[I|)). As voltage is applied, electrons incoming from the left with energies such that e nL /e < V will 
see a potential step-up. The number of these states is given approximately as 



V2em*V/hAk. (16) 



For states incoming from the left, we approximate T; ~ for incoming energies less than the potential step height, 
and Ti ~ 1 for energies greater than the potential step height. In contrast, electrons incoming from the right see a 
step-down potential and we approximate T r ~ 1 for all electrons incoming from the right. The electron current can 
then be estimated as 
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m*L 

n 

n=l 
ny 



m*L 

n=l 
, nv 

eh 



m*L 

n=l 

with the convention that current flow is opposite the direction of electron flow. For small Ak and large np (these are 
standard conditions for derivation of the Landauer formula) , we have 

eh 



Ak ndn=-V. (18) 
m*L J h 

The current and voltage yield a conductance go = e 2 /h and the Landauer result for conductance quantization is 
obtained. 

This derivation seems odd, as it appears that it is not the states at the Fermi level that contribute to the current, 
but states low in energy (or momentum) that yield a current. The situation can be summarized in fig. E]. In fig.[4^), 
the product of the momentum occupation and the transmission coefficients for left and right states are given, whereas 
in fig. |4Jd), the corresponding product of the energy level occupations and transmission coefficients are given. In 
fig. Bk), it appears the currents arise from low momentum states, in fig. Bp), currents appears to arise from energy 
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states at ep. However, there is no contradiction. If the currents from momentum states with energies e < ep are 
summed, the currents associated with states below ep cancel. 

The current incoming from the left can be rewritten as an integral over energy 

e 2 f ep 

h = — r / de, (19) 

n JeV 

using e = and de = -^kdk. Similarly the current from the right is re- written but now with e = + eV 



e 2 r^F+eV 



r ~ h 

The resulting current is 



de. (20) 

eV 



2 rtF+eV 2 

'si < 21 > 

which will be recognized as the more familiar form for expressing conductance quantization. The physics is consistent 
whether calculating currents using the momentum distributions fig. Ea) or the energy distributions fig. Er 



V. DISCUSSION 



The first point of our presentation is that the MECS approach is consistent with a Landauer description of electron 
transport. BK's conclusion ascribes the failure of their calculations for a single particle and a correlated model to 
the MECS formulation of boundary conditions. The analysis of sect. |III| highlights this is not the case; based on an 
analytical model, the use of the MECS formulation is consistent with a scattering approach to electron transport. 
The advantage of an analysis based upon an analytical model is that it avoids issues associated with linear response 
approximations, perturbation theory, variational methods in a finite basis, specific implementations of electronic 
structure, or other numerical approximations, thereby allowing a clear focus on the physical assumptions made when 
using the method. In sect. |IV[ we continue in this vein and demonstrate how the MECS formulation reproduces 
conductance quantization in a system of non-interacting electrons scattering off a potential barrier. 

We would also like to touch upon some formal points raised in the BK work. The authors criticize the use of a 
configuration expansion to describe transport problems. We note that it is shown in several works that a properly 
designed variational calculation can provide accurate properties governing electron transport such as electronic spec- 
tra [TB] with compact expansion vectors [T7]- The variational structure of MECS calculations performed to date has 
also been noted by BK and we believe misinterpreted. Using a variational approximation to the wavefunction does 
not result in an exact eigenfunction of the system Hamiltonian but rather the best approximation using the functional 
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< \E'|if| , I' > / < V E'| V E' > for the approximating function \*f? > and subject to the application of the constraint condi- 
tions. As a consequence, it is well known that integrated quantities such as the energy may be better approximated 
compared to local properties such as spin density or electron current density. Hence in previous MECS calculations 
the possibility to introduce constraints to enforce current conservation on a scattering region was introduced. How- 
ever, this is a numerical feature related to the variational nature of the calculations and the application of the open 
system boundary conditions does not imply the violation of current conservation contrary to a supposition in BK. 
We have already noted the origin of the current variations from variational calculations. It also well known that, for 
example, perturbation theories do not conserve many physically conserved quantities, including electronic current. 
Considerable care is needed in defining finite expansions that are current conserving |18j . Hence we maintain the 
BK have incorrectly ascribed to MECS the curent oscillations they calculate in a tight binding model within linear 
response to the boundary conditions applied within MECS. 

As a final point, the application of the Wigner constraints in the case of a one-dimensional problem as we have 
introduced and as attempted for a tight-binding linear chain in BK requires particular care in the following sense. 
Formally, for a single electron model of metallic electrodes (a reasonable assumption), the density matrix decays as 
l/\x — x'\ d with d the spatial dimension of the electrode. How one decides to treat this long range behaviour influences 
the boundary conditions. Another way to express this is that the density matrix does not satisfy Kohn's principle 
of 'nearsightedness' for these examples [19], whereas in three-dimensional models of metal electrodes [TJ [3j [lOj [11] 
the density matrix decays within typically less than 0.5 nm |20] thereby greatly simplifying the introduction of open 
system boundary conditions through use of the Wigner function. The decay of the density matrix in three-dimensional 
metallic electrodes helps in the calculation of the equilibrium (V—0) Wigner function with small explicit electrode 
regions, without coupling between the electrodes or coupling the electrode equilibrium regions to the scattering region. 

VI. CONCLUSION 

We have shown that the MECS boundary conditions and introduction of a chemical potential imbalance between 
electrodes reduces to the correct single particle limit, and hence the claim in BK for the failure of the MECS boundary 
conditions to reproduce conductance quantization is incorrect. The model analysis provided reveals that a failure of a 
calculation attempting to apply the boundary conditions using the Wigner function for momentum distributions for 
this or related models can not be attributed to a failure of the MECS formulation. 

We have also provided an analysis on both the formulation of many-electron scattering using the Wigner function 
boundary conditions, and touched upon issues related to the numerical implementation of the model. We will present 
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a similar analysis on more realistic models of atomic and molecular scale systems and consider the effect of various 
numerical approximations on MECS transport calculations in future work. 
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FIG. 1: The Fermi-Dirac energy dispersions for incoming electrons from two electrodes described by parabolic bands, a) with 
no voltage difference between the electrodes, b) with applied voltage, c) The corresponding momentum distribution functions 
for a) and b). Note that there is no difference to the momentum distributions with application of voltage. 
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FIG. 2: A simple model of the action of two electrodes in generating an electric field. The application of a chemical potential 
difference results in surface charges generating an electric field between the electrodes. 
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FIG. 3: The Wigner function calculated with scattering wavefunctions defined on a region [-L/2, +L/2] and with normalization 
chosen such that a completely occupied state is with / = 1. a) The Wigner function calculated at qi = — L/4 for a region of 
length L — 2 nm with /cf=0.12 nm _1 =0.635 a.u (corresponding to gold electrodes) with V= (blue), 1 (purple), 2 (green) 
volts, b) Same as in a) but with L=20 nm. c) The Wigner function calculated at q r — +L/4, otherwise the same as in a), d) 
Same as in c), but with L—20 nm. 
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FIG. 4: a) The momentum distributions corresponding to the model of sect. IV b) The energy distributions for the model of 
sect. |IV| Both distributions may be used to calculated conductance quantization. 



